Accelerator dynamics of a fractional kicked rotor 
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It is shown that the Weyl fractional derivative can quantize an open system. A fractional kicked 
rotor is studied in the framework of the fractional Schrodinger equation. The system is described 
by the non-Hermitian Hamiltonian by virtue of the Weyl fractional derivative. Violation of space 
symmetry leads to acceleration of the orbital momentum. Quantum localization saturates this 
acceleration, such that the average value of the orbital momentum can be a direct current and the 
system behaves like a ratchet. The classical counterpart is a nonlinear kicked rotor with absorbing 
boundary conditions. 
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Application of fractional calculus to quantum processes 
is a new approach to the study of fractional properties 
of quantum phenomena 0, 0, IS S IE 01 ■ In this Report 
we consider a quantum chaotic dynamics of a fractional 
kicked rotor (FKR) . The Hamiltonian of the system is 



H = T + ecosx S{t — n) , 



(1) 



where e is an amplitude of the periodic perturbation 
which is a train of 6 kicks. The kinetic part of the Hamil- 
tonian is modeled by the fractional Weyl derivative 



f ={-ih)°'W°'/a., 



(2) 



where /i is a dimensionless Planck constant, and a = 2-/3 
with < ;9 < 1. When a = 2 Eq. (HJ is the quan- 
tum kicked rotor J^. For a periodic function f{x) = 
^/fcC^**"'^, the Fourier transform property determines 
the fractional Weyl derivative W" in the following sim- 
plest way (see |3j, ch. 4.3) 



W"/(x) 



(-zA:)"/fee- 



(3) 



Since only periodic functions are considered here, this 
oversimplified definition is sufficient without the bur- 
den of fractional calculus details 8]. Thus, the kinetic 
term in the Hamiltonian is defined on the basis 



|fc) = e'''^/v27r as follows 



f\k)^T{k)\k) 



(hk) 



2-/3 



2-/3 



-\k) 



(4) 



This non-Hermitian operator has complex eigenvalues for 
fc < 0, which are defined on the complex plain with a 
cut from to —oo, such that 1^^ = 1 and (—1)"'' = 
cos /Jtt — i sin /Jtt, and therefore, k^^ — \k\^/^e^"^^''\ 
where /3(fc) = (3[1 - sgn(fc)]/2 It is worth men- 

tioning that the fractional derivative in Eq. ^ appears 
naturally in quantum lattice dynamics with long range 
interaction Q, where {—ik)" is a particular case of poly- 
logarithm (see Appendix in Ref. p|). 



A quantum map for the wave function tl>{x,t) is 

'4!{x,t + l)^U%l>{x,t), (5) 
where the evolution operator on the period 



tj = exp 



-le cos a; 



exp 



-iT 
h 



(6) 



describes free dissipative motion and then a kick. Dy- 
namics of the FKR is studied numerically, where Eq. 
enables one to use the fast Fourier transform as an ef- 
ficient way to iterate the quantum map lO). A specific 
property of this Hamiltonian dynamics is quantum dissi- 
pation resulting in probability leakage and described by 
the survival probability 



(7) 



where |/np is probability of level occupation at time t. 
The initial occupation is fn{t = 0) = 5n,Q- Another spe- 
cific characteristic is the nonzero mean value of the or- 
bital momentum 



(8) 



due to the asymmetry of the quantum kinetic term T . 
Results of the numerical study of the quantum map © 
are shown in Figs. 1-4. The quantum dissipation leads 
to the asymmetrical distribution of the level occupation 
\fn[t)\ (see Fig. 1) that results in the nonzero first mo- 
ment of the orbital momentum (p) ^ V'^ in Fig. 2. Quan- 
tum localization saturates the acceleration with time. 
This accelerator dynamics is accompanied by the power 
law decay of the survival probability P{t) ^ t~'*^ with the 
exponent 72 ~ 0.71 shown in Fig. 3, and then the de- 
cay rate increases with the time due to quantum effects. 
Quantum localization affects strongly both 71 and 72 . By 
increase of the quantum parameter, when h — 0.76, the 
exponent 71 approaches zero (in Fig. 4 the slope is 10~^), 
and the survival probability decays at the rate 72 « 0.99. 
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FIG. 2: Acceleration of the average orbital momentum for the 
same parameters as in Fig. 1. The insert is the log- log plot, 
and the solid line corresponds to 71 — 0.35 obtained by the 
least squared calculation. 



FIG. 4: Quantum saturation of {p{t)) due to localization when 
h = 0.76, P = 0.05 and the same e as in Figs 1-3. The slope of 
the solid line is 10~^. The insert shows the power low decay 
of P{t) with 72 ~ 0.99 due to the linear interpolation. 



To understand the obtained numerical results and the 
physical relevance of the fractional Schrodinger equation, 
the classical limit h ^ is performed in the Wigner rep- 
resentation. Thus, the system is described by the Wigner 
function W{x,p, t) which is a c-number projection of the 
density matrix in the Weyl rule of association between 
c-numbers and operators. The Weyl transformation of 
an arbitrary operator function G{x,p) is [9I flo| 



F{x,p) =Tt [G{x,p)A{x -x,p-p)] , 



(9) 



where F{x,p) is a c-number function and A(x—x,p—p) is 
a projection operator which acts as the two dimensional 
Fourier transform. For the cylindrical phase space the 



projection operator is pT| 

00 ^ 



A(x- 



im{x — x)-\-i^ (p—p) 



(10) 

This operator determines an inverse transform as well: 



Gix,p) 



where p - 
pit) is 



00 

E 

k— — oc 



— I F{x,hk)A{x^ 



x,hk — p) , (11) 



hk. The quantum map for the density matrix 



p{t + l) = U^p{t)lJ . 



(12) 
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Therefore, evolution of the Wigner function 

W{t,x,p)= Tr [p{t)Mx ~x,p~p)], 
for the period determined by the map H12(l . is 

W(t + l,x,p) ^Tt U''p{t)UA{x-x,p~p) ] 
= E^=-ooIo^ Kj^i^,PW,P'Wit,x',p')dx' , (13) 
where Kj:^(x,p\x' ,p') is Green's function for the period 
Kf^{x,p\x',p') = 2^ e™(---'+«')e'«'(p'-rt 



exp ^T*{p + hm/2)- ^T{p-hm/2) 



X exp 



law phenomenon due to quantum tunneling has been an 
issue of extensive studies in quantum chaos • 

Quantum localization leads to the exponential restric- 
tion of the initial profile spreading in the orbital mo- 
mentum space from above. This property results in sat- 
uration of acceleration of {p{t)); namely, at t — > oo it 
follows that {p{t)) const. Such a behavior is found for 
h = 0.76 and /3 = 0.05. In Fig. 4 one sees a direct current 
of {p{t)) for 5 • 10^ iterations and K = 3. This double im- 
pact of asymmetric absorption and quantum localization 
leads, asymptotically, to a quantum like ratchet which 
differs from quantum one obtained on a classical chaotic 
attractor [l^ . 

It is worth mentioning that in the class of periodic 
functions, eigenvalues of the unperturbed Hamiltonian 



f cos{x' + hC/2) ~ I cos(x' - /i^72)J . (14) r coincide with Fo(p) = (-«^|;)": and have the same 



The trace is Tr[...] = EfcC*'!---!^)- '^'^ classical limit 
/i ^ 0, we obtain in Eq. 1)14(1 that the difference of the 
perturbations in the exponential is — iecosx, while the 



difference of the kinetic terms is imp 



1-/3 = 



){p) for 



p > and -2sin(/37r)T(|p|)//i for p < 0. The last term 
diverges at /i = and yields identical zero for Green's 
function Kj^^^j{p < 0) = 0. Thus, the classical Green 
function 

^h=oi^iPWp') — Q{p)S{x — X — uj(j)))5{p — p— esina:') 

(15) 

corresponds to the classical map M. 

Pn+l = Pn + esina;„, Xn+l Xn + (^(Pn+i) (16) 

of the nonlinear kicked rotor with the nonlinear frequency 
ijj{p), and absorbing boundary conditions for p < 0, that 
the Heaviside function Q{p) reflects. 

Therefore, the fractional Hamiltonian Q corresponds 
to the open system Eqs. (|15|) . l|16|l . Chaotic dynamics of 
this open system takes place in the upper half of the cylin- 
drical phase space. The stability property is determined 
by the trace of the linearized map dM 



Tr [dT] = 2 + e(l - f3)p'^ cos x . 



(17) 



For any e there are stable regions {Ax, Ap} determined 
by the locus of elliptic points {xe,Pe} (see Fig. 5) 



Xe = arccos 



2pf 



e(l-/3) 



(18) 



The presence of this infinite regular elliptic island struc- 
ture, which leads to the stickiness of chaotic trajectories, 
also results in the power law decay of the survival proba- 
bility for the quantum counterpart in Fig. 3. This power 



classical limit of Eq. H15|l . This local derivative has the 
classical counterpart with the Hamiltonian Hq{p) = p" 
which does not coincide with Eq. p5|l . Namely, the 
Hamiltonian Ho{p) is the classical system with dissipa- 
tion for p < 0, while the map M in Eqs. H15|l and H16() 
is the open system where a particle is set apart from the 
dynamics for p < 0. 

Fractional Schrodinger equation 



(19) 



describes quantum dissipative Hamiltonian dynamics. 
The classical counterpart is a nonlinear motion with dis- 
persion uj{p) realized on the upper half plain of the phase 
space with absorption in the lower half plain. It has 
well defined physical meaning. Therefore, the fractional 
Schrodinger equation 1(19(1 can be a generalized approach 
for any functions for which the Fourier transform is valid. 
In this case, the opposite classical~to-quantum transition 
can be performed by determining the Heaviside function 
in Eq. lO 



Q{p) = lim exp 



where T(p) is uniquely defined by the condition w(p) = 
T'{p). Thus, fractional derivatives quantize classical 
open systems in the framework of the non-Hermitian 
Hamiltonians. 
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